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Abstract: 
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factors observe several relations inspired by heavy-quark symmetry. 
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1 Introduction 



The challenge to understand the physics of CP- violation related to the structure of the CKM 
mixing matrix in (and beyond) the Standard Model is fuelling an impressive experimental 
programme for the study of S-decays, both exclusive and inclusive. Abundant data in various 
exclusive channels are expected to arrive within the next few years from the dedicated B- 
factories BaBar and Belle; their potential impact on our understanding of CP-violation at 
the electroweak scale will crucially depend on our possibility to control the effects of strong 
interactions. For exclusive decays with only one hadron in the final state, the task is to 
calculate various transition form factors; it has already attracted significant attention in the 
literature. 

In this paper we present the first complete results for the exclusive semileptonic and rare 
radiative 5-dccays to light vector-mesons in the light-cone sum-rule approach. Exclusive 
decays, which are the principal concern of this work, can be grouped as semileptonic decays: 

• B^,d pev, 

rare decays corresponding to 6 — > s transitions, which we term CKM-allowed: 
. Bu,d ^K* + -f, Bu4 ^K* + Z+r , 

• 5, ^ + 7, Bs ^(f) + l+l-, 

and b ^ d transitions, which we call CKM-suppressed: 

• -Bd ^ (p, a;) 7, B^ ^ {p,uj) + l+l' , 

• -B« ^ P + 7, p + 1+1-, 

• Bs^K* + -f, Bs^K* + l+l-. 

Let y be a vector-meson, i.e. p, ur, K* or (p, and let p^, e* and mv be its momentum, 
polarization vector and mass, respectively. Let ps {rriB) be the momentum (mass) of the 
S-meson. We define semileptonic form factors by {q — Pb — p) 



{V{p)\{y -A)^\B{pb)) = -iel{mB + mv)A\{q')+i{pB+p)^,{,e*pB 



niB + my 



+ Me*PB) ^ {Aliq^) - <(g^)) + 6,.,.6*>^p'^ . (LI) 



Note the exact relations 



y 2\ mB+my .y, 2\ nis — rXly .y, i\ 

<(o) = ^r(o), 

{V\d^A'^\B) = 2my{e*pB)AUq'). (L2) 
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The second relation in ( p. .21 ) ensures that there is no kinematical singularity in the matrix 
element at = 0. 

Rare decays are described by the above semileptonic form factors and the following 
penguin form factors: 

{V\i:a^,q''{l + j5)b\B{pB)) = ze^,,.e*Xp^ 2ri(g2) 

+ T^iq') {e;(m| - m^) - (e^^) {ps + p),} 

+ T^{q^){e*pB) [q^ - , , {Pb+P),\ (1-3) 



2 2 

— rriy 



with 



Ti(0)=T2(0). (1.4) 



Here ip = s,d. All signs are defined in such a way as to render the form factors positive. 

The physical range of extends from g^j^ = to q'^^^ = {niB — mvY for three-body 
decays and = for two-body decays. 

The method of light-cone sum-rules was first suggested for the study of weak baryon- 
decays in and later extended to heavy- meson decays in [Q]. It is a non-perturbative 
approach, which combines ideas of QCD sum-rules ||^ with the twist-expansion characteristic 
of hard exclusive processes in QCD [Q and makes explicit use of the large energy of the final- 
state vector-meson at small values of the momentum-transfer to leptons, g^. In this respect, 
the light-cone sum- rule approach is complementary to lattice calculations [Q, which are 
mainly restricted to form factors at small recoil (large values of g^). Of course, the light- 
cone sum-rules lack the rigour of the lattice approach. Nevertheless, they prove to provide a 
powerful non-perturbative model, which is explicitly consistent with perturbative QCD and 
the heavy-quark limit. 

Early studies of exclusive 5-decays in the light-cone sum-rule approach were restricted to 
contributions of leading-twist and did not take radiative corrections into account, see Refs. P, 
1^ for a review and references to original publications. Very recently, these corrections have 
been calculated for the semileptonic B n^Kev decays In this work we calculate 
radiative and higher-twist corrections to all form factors involving vector-mesons (see above) 
making use of new results on distribution amplitudes of vector- mesons, reported in |10|, |ll|]. 
We find that the corrections in question are fairly small in all cases. 

The presentation is organized as follows: in Sec. 2 we recall basic ideas of the light-cone 
sum-rule approach and derive radiative and higher-twist corrections to the form factors in 
question in a compact form. Section 3 presents our main results and includes a discussion of 
input-parameters as well as error estimates. In Sec. 4 we discuss relations between semilep- 
tonic and penguin form factors in the heavy-quark limit. Section 5 is reserved to a summary 
and conclusions. The paper has two appendices: in App. A we collect the relevant loop 
integrals for the calculation of radiative corrections. Appendix B contains a summary of the 
results of |10|, |lT| on vector-meson distribution amplitudes. 
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2 Method and Calculation 



2.1 General Framework 

Consider semileptonic peu and rare — > K*i~^i~ decays as representative examples. 

We choose a i?-meson "interpolating current" js = dij^b, so that 

{0\jb\B{pb)) = ^, (2.1) 
rrib 

where fs is the usual -B-decay constant and the 6-quark mass. In order to obtain 
information on the form factors, we study the set of suitable correlation functions:^ 



+ tV^pl, q') ^ (g + 2p), + zr-(p|, q') ^ + r^(p|, q^)e;^P^ e^q^p,, (2.2) 
I I d'ye-'P-y{K*{p)\T[-sa,,^,bmjUy)\0) = A{pl,q'){e;{2p + q), ~ e:{2p + q),} 

- I3{pl, q^){elqy - e*g^)} - 2C(p|, q^) — {p^q^ - q^p^}. (2.3) 

The Lorentz-invariant functions F*^'^'^, A, B, C can be calculated in QCD for large Euclidian 
p\. More precisely, if —Ps'^^i then the correlation functions (|2.2| ), (p.3|) are dominated 
by the region of small y"^ and can be systematically expanded in powers of the deviation from 
the light-cone y^ = 0. The light-cone expansion presents a modification of the usual Wilson 
operator product expansion, such that relevant operators are non-local and are classified 
in terms of twist rather than dimension. Matrix elements of non-local light-cone operators 
between the vacuum and the vector-meson state define meson distribution amplitudes [^], 
which describe the partition of the meson-momentum between the constituents in the infi- 
nite momentum frame. In particular, there exist two leading-twist distribution amplitudes 
for vector-mesons, see App. B, corresponding to longitudinal and transverse polarizations, 
respectively: 

{p\u{0)^^d{z)\0) = fpm.p/— fdue'^P^<P\\{u,p), (2.4) 

pz Jo 

(p|fZ(0)a^,rf(z)|0) = -^fJ(p)(e;p,-p^e:) / rfw e^"^^ 0x(w, /i), (2.5) 

Jo 



and similarly for K* and 0. Here z is an auxiliary light-like vector, u is the momentum 

fp^ fp 



fraction carried by the valence quark, and the decay constants fp, fj are defined in App. B, 



^ In this work we define invariant functions with respect to the Lorentz-structure ^ instead of e*(j |1S 
in order to remove a kinematical singularity for p 0. 
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while /i specifies the scale: extraction of the leading asymptotic behaviour in field theories 
invariably produces singularities, which reflect themselves in the scale dependence of dis- 
tribution amplitudes. As always, this scale dependence cancels in physical quantities by a 
corresponding dependence of coefficient functions. 

The invariant amplitudes in ( |2.2D , ( ^.31 ) can be calculated in terms of meson distribution 
amplitudes, in complete analogy with the calculation of structure functions in deep inelastic 
lepton-nucleon scattering in terms of nucleon parton-distributions: the off-shellness ml — p\ 
plays the role of photon virtuality Q^. As an illustration, consider the tree- level leading-twist 
result for r°, adapted from Ref. ||T2| : 



Jo mi - ^Pb - 



. -,.Udu)'\^. (2.6) 

We want to emphasize that the procedure is rigorous at this point: all corrections can (in 
principle) be included in a systematic way, and their evaluation is precisely the subject of 
this work. 

The subtle part concerns the extraction of the S-meson contribution to the invariant 
amplitudes. The exact amplitude r° (in nature) has a pole at p% = m\ corresponding to 
the intermediate 5-meson state, and this contribution can be written in terms of the form 
factor A^~*'' defined in (|1.1[): 

r°B-meson = {^b + m,)A^r''{q') " " ^ • (2-7) 

On the other hand, the QCD calculation at <^ ml is only approximate and, continued 
analytically to "Minkowskian" > ml, produces a smooth imaginary part with no sign 
of a pole behaviour. To proceed, we invoke the concept of duality, assuming that the exact 
spectral density and the one calculated in QCD coincide on the average, that is integrated 
over a sufficient region of energy. In particular, we assume that the 5-meson contribution 
is obtained by the integral of the QCD spectral density over the duality region: 

so 

\ f d 

r^-meson = ^ / Disc,|r°QCD(^, (2-8) 

ml 

The parameter sq ~ (34-35) GeV^ is called "continuum-threshold" and is fixed from QCD 
sum-rules for fs, see e.g. ||13[. Equating the above two representations, one obtains a light- 
cone sum-rule for the form factor Ai. Sum-rules for the other form factors are constructed 
in precisely the same manner. 

While the accuracy of the QCD calculation can be controlled (and improved), the duality 
approximation introduces an irreducible uncertainty in predictions for the form factors, which 
is usually believed to be of order (10-15)%. Practical calculations in the sum-rule framework 
involve some technical tricks to reduce this uncertainty, e.g. Borel transformation, which we 
will not discuss here. These techniques are well established and their detailed description 
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a b 






Figure 1: The leading-order diagram (a) and one- loop radiative corrections (b-f). 



in the particular context of light-cone sum- rules can be found, for instance, in Refs. [0, |T^] . 
The work [|T2| also contains a detailed comparison of the light-cone sum-rule approach to 
traditional QCD sum-rules and can serve as introduction for the more theoretically-minded 
reader. 



2.2 Radiative Corrections 

Radiative corrections to the sum-rules correspond to one-loop corrections to the coefficient 
functions in front of leading-twist distribution amplitudes and are given by the diagrams 
shown in Fig. 1. The calculation is done in dimensional regularization, and it is sufficient 
to consider matrix elements over on-shell massless quark and antiquark carrying momentum 
fraction up and up, respectively. The transversely polarized and longitudinally polarized 
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meson states are projected on by 



{V±{p)\ua{0)di,{x)\0) = -l/^[a^,],„e>-y due'^'^'^Mu) (2-9) 
= -I fv[<y^''l^]bae^.pae*''p'' j du e'^P>x(n), (2.10) 
{V\\{p)\uMdb{xm = - fvmv[T/>],a— j rfne™P^0,|H 

4J: 'PX J 

m2 1 f 

^ -fvWa due''^^-ct^\\{u), (2.11) 



where a, h are spinor indices. In the last hne in (|2.11|) we made use of the fact that for ultra- 
relativistic longitudinal vector-mesons Viil°^v up to 0(my/|p|^) corrections. This is a 
justified approximation for the calculation of radiative corrections to leading-twist accuracy, 
to which end the meson- mass can be neglected throughout. For further use we introduce the 
notation for the projection operators: 

P± = -\ilo^pe*^pP or Vf = -l/^a„;375e"^^^e>.. (2.12) 

In what follows, they will be treated as D-dimensional objects . 

The calculation in question is in principle straightforward and similar to the existing 
calculations of NLO corrections to hard exclusive processes |]T^. One has to consider one- 



loop diagrams with a heavy-quark and two different kinematic invariants and p%, which 
makes formulas rather cumbersome. The specific requirement is to organize the expressions 
in a form suitable for a dispersion representation in p^, cf. Eq. ( p.8|) , so that continuum 
subtraction can be made. 

Analytic expressions for 5-decays to light pseudoscalar-mesons vr, K have been made 
available recently For vector-mesons the number of form factors is so large that working 
out (relatively) compact analytic expressions is not worth the effort. In this work we prefer to 
give the formulae in terms of traces and general momentum integrals (see below and App. A), 
which can be compiled and evaluated numerically using the mathematica programming 
language^. 

A usual subtlety concerns the treatment of 75. The results for the form factors given 
below are obtained using "naive dimensional regularization" (NDR), and the same scheme 
has to be applied to the calculation of Wilson coefficients for penguin operators. 

There are two form factors in whose calculation one encounters an odd number of 75 
in traces, which could cause ambiguities: V and Ti. Only transverse mesons contribute 

^ The computer code is available from P.B. upon request. 
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to these form factors. In both possible ambiguity comes solely from the 5-vertex 

correction in Fig. Id, whereas in all other diagrams contraction of 7 matrices over 75 can 
be avoided. There are several ways out: (a) use a 't Hooft-Veltman prescription for 75 and 
apply a finite renormalization to restore the Ward identities, as in |T^; (b) instead of the 
"natural" projection ( |2.9D , use ( ^.101 ), which introduces a second 75 and thus eliminates the 
problem; (c) modify the definition of the form factors ( |1.3| ) to 

{V\m^,^,h\B) = A{q^) {e;(pB - (p^ - B{q^) {t*^q, - q^tl) 



-C{q' 



m 



B 



7717 



{(Pb + p)^,qu - q^{pB + p)u} ■ 



Using 



a 



pa 



and contracting with q", one finds 

A{q')=T,{q% 



(2.13) 



771 



B 



Uq') - 



[Tr{q')~Uq' 



mt 



[T,{q')-Uq')] 



(2.14) 



from which the relation ( |1.4|) follows. It is thus sufficient to calculate A, B and C instead 
of Ti with the premium to avoid any 75 problem. We have checked that all of the above 
prescriptions yield identical results. 

After these preliminary remarks, we are now in a position to calculate the diagrams in 
Fig. 1. The tree- level contribution of Fig. la equals 



rpio) ^ 1 Tr(r(^5 -u^ + mb)75P), 



(2.15) 



where F is the Dirac structure of the weak vertex, V is one of the projection operators defined 
in Eqs. (|2l^ ), and 



— uq'^. 



It proves convenient to replace in Eq. ( p.l5| ) the running MS 6-quark mass by the one-loop 
pole mass, which is given by 



771 



pole 



47r2 



This replacement induces the radiative correction 
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47r2 



1 



3 , 771 

In —I 



4 



3 

- In — 

4 /i2 



Tr('PF(^^ - + mfe)75 



(2.16) 



- mfeCF — 
s 47r^ 



3 rn^ 

1 - 4 In -r ) MVWb 



u 



(2.17) 
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The general strategy is to simplify the traces as much as possible, but to keep V and V 
arbitrary. Also contraction of 7 matrices over 75 is only allowed in the 5-vertex correction. 
It turns out that all one-loop diagrams can be expressed in terms of the following traces: 





= 'iv(prg^75) = 


1^ (^rj^BTs) , 




= Tr(Pr75), 






= IV(7'^r75) = 


T^Wsr75), 


IV4 


= Tr(P^r^75). 





(2.18) 



Let us also introduce 

arV := 7.^7", arT := 7,17". (2.19) 
The 6-quark self-energy diagram in Fig. lb is: 

ySE = { [4.ml {Y + {1- e)Z) + 2s{l - e){Y - Z)] ( Tr^ + m.Tr^) 

-2mfes(r + (l-e)Z) Tra}, (2.20) 

where D = 4 — 2e and the expressions for momentum integrals Y, Z are given in App. A. 
The self-energy insertions in external light-quark legs in Fig. Ic only contribute logarithmic 
terms in dimensional regularization: 

ySE ^ ^ + ^ _ In , (2.21) 

where we distinguish between the ultra-violet scale Huv, which is to be identified with 
the renormalization scale of the curent Jb and the penguin operators, and the infra-red 
renormalization-scale /im corresponding to the factorization scale in meson distribution am- 
plitudes. 

For the i?-vertex correction in Fig. Id, one obtains: 

T^^2 { (-8(7(1 - e) - 1 - mlB + u{pl - q^)A) ( Tri + m,, Tr^) - m^sB Tra} . 

^ (2.22) 
For the weak vertex in Fig. le we find: 

rpw ^ _ y__F_ c ( Tn + mft Tr2) + av D {q" Trg + mb Tr4) + {pi - q^)uar E Trg 
s 

+2(p| - q')uA{ Til + rrib Tr^) + m^B [2(m6 IVi + TVs) - ar(m6 IVg + IV4)] } . (2.23) 

Finally, the box-diagram in Fig. If can be written as 

ybox _ ^g-^CFar {arH{ Tn + m, T12) + /(-m^ Tr4 + {s - ml) Trg) + Trg) , (2.24) 

where -B„=i is the limiting value of B for m — > 1. 

Definitions and explicit expressions for the one-loop integrals A, B, C, etc., are given in 
App. A. 
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a 



b 



Figure 2: The higher-twist contributions. 



2.3 Higher-Twist Contributions 

Higher-twist terms generically refer to contributions to the hght-cone expansion of the cor- 
relation functions ( |2.2| ) and (P73|), which are suppressed by powers of — p^). In the 
sum-rules, such corrections are suppressed by powers of the Borel parameter. Higher-twist 
corrections are usually divided into "kinematical" , originating from the non-zero mass of the 
vector-meson, and "dynamical" , related to contributions of higher Fock-states and transverse 
quark-motion. In this paper we take into account both effects to twist-4 accuracy, making 
use of the new results on distribution amplitudes of vector-mesons reported in 
summarized in App. B. 



10, 111 and 



The calculation is most conveniently done using the background-field approach of [|1^ 
The diagrams of the type shown in Fig. 2a are taken into account within this method by the 
expansion of the non-local quark-antiquark operator in powers of the deviation from the light- 
cone; they give rise to contributions of two-particle distribution amplitudes of higher- twist, 
see Eqs. ( B.12| ) and ( B.27 ). The contribution of the gluon-emission from the heavy-quark is 



calculated using the light-cone expansion of the quark-propagator [jT6|, p!7| : 



{0\T{b{x)m}\0) 



iSb{x) - ig 



(2^ 



-ikx 



dv 



1 It + TTlb 

2 {ml - P)2 



m, 



,(2.25) 



where Sb{x) is the free quark-propagator. As in the case of radiative corrections, our strategy 
in this work is to derive the most general expression for all form factors in question, suitable 
for implementation in analytic/numerical calculations using mathematica. We obtain 



CF = - / du \ ifyfriv 



II dQo, pq 16 dQpdQp 



Tr{TSb{Q)i, 



pq 2 dQa 



9 



-7 e«/57^e*V^7fV) T^Tr(r5,(Q)7o 
4 (7(^5 



T 



1 (9^ 

myAxiu) 
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d_ 

dQc 



dQfs 



Tr{TSb{QHcr. 



Q/3j 



+-/vmy y dv jVa 



2v{pq) {A{a) + V{a))Tr{Tf 



+ ^ ( 2$(a) + ^(a) - 2<l>(a) - ^(a) ) Tr(r(^ + mfe)^75j 
pq 



+4m^yv{e*q) $(a) - $(a) Tr(r75) - m^^t; — ^(«)Tr(r^^7 



+-fymy I dv jVa ml - {q + {ai + va3)pf 



2v{t*q)T{a)Tr{Ti>^^ 



+ (5(a) - S{a) + - T^la) + T^t^a) - T^' {a) ) Tr{T{4 + m,)fi>^,) 



.(4) 



(4)/ 



.(4), 



+ 2v (rf (a) - Tt\a) - 5(a) - 5(a)) [(e*g)Tr (1^75) - (pg)Tr (r/*75) 



+ 2t;(Ti')(a)-Tf^(a)-5(a 



(4)/ 



75 J 



(2.26) 



where Q = q + up and CF G {r°'''='^, i3, C}. Definitions and explicit expressions for tlie 
numerous distribution amplitudes are collected in App. B Q In addition, we use the notation 



$PfM) 



(2.27) 



To leading-twist accuracy, our result agrees with the expressions available in the literature, 
18|, 0, 



see 



■^Despite appearance, the number of non-perturbative parameters in the description of higher-twist dis- 
tributions is small, since they are related by exact equations of motion, see jl^, |l^ and App. B. 

^ The sum-rule for Ti given in |l^ misses a contribution of $|| ; this term can be formally viewed as 
part of the kinematic higher-twist correction, which is included in |l^ only partially. 
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3 Results 



In this section we present results of the numerical analysis of the light-cone sum-rules for 
the form factors defined in ( |1.1| ) and ( p. .31 ) for B- and -B^-decays. The sum-rules depend 
on several parameters, those needed to describe the 5-meson on the one hand and those 
describing the vector-meson on the other hand. The former are essentially J'b (/bJ, the 
leptonic decay constant defined in ( |2.1| ), the 6-quark mass rrih, the continuum-threshold Sq 
introduced in ( p.8D , and the Borel parameter mentioned in Sec. 2.1. Lacking experimental 
determination of Jb and Jb^, we determine their values from two-point QCD sum-rules to 
0{as) accuracy (see e.g. |jl3|); this fixes sq, which depends on nib, and also the "window" in 



in which the sum-rules are evaluated. We then use the same values for rrib, sq and 
in both the QCD sum-rule for J'b and the light-cone sum-rules for the form factors]^ which 
helps to reduce the systematic uncertainty of the approach. The corresponding parameter- 
sets and results for the decay constants are given in Table |l|. The question of the value of 
the 6-quark mass has attracted considerable attention recently; following these developments 
20[] , we use the value rrib = (4.8 ± 0.1) GeV. Our results for J'b agree well with new lattice 



determinations [pi 



The parameters of light mesons are collected in App. B, Tables |^ and 0. These parame- 
ters are evaluated at the factorization scale fijj^ = m?^ — "^l ~ 4-8 GeV^, which is the typical 
virtuality of the virtual 6-quark in the process. The penguin form factors also depend on 
the ultra-violet renormalization scale of the effective weak Hamiltonian, for which we choose 
Huv = ^b- Using the central values of all parameters, we obtain the form factors plotted 
in Figs. 1^ and ^. For their representation in algebraic form, a parametrization in terms of 
three parameters proves convenient: 

Fi<l') = . ..2 . (3.1) 



l-ap^ + bpl- 



with the fit parameters F{0), ap and bp- Here rriB is the mass of the relevant 5-meson, i.e. 
TfiB^^ for i?„ rf-decays and for S^-decays. This parametrization describes all 28 form 
factors to an accuracy of 1.8% or better for < < 17GeV^. 

Let us now discuss the dependence of the results on the input-parameters and approxi- 
mations involved. First we note that the net impact of radiative corrections is very small at 
small and at most 5% at = 0. Their effect increases, however, at large and leads to 
a decrease of the form factors A2 and T3 at = 17GeV^ by 20% with respect to their tree- 
level values; the impact on the other form factors stays in the 5% range. The small effect of 
radiative corrections was anticipated in the tree-level analysis of Ref. [jl2| and also observed 



in the calculation of 0(0;^) corrections for B — > pseudoscalar decays 0. It is due to the 
fact that the biggest contribution to radiative corrections (in Feynman gauge) comes from 
the S-vertex correction diagram, which enters both the calculation of Jb and the light-cone 



^To be precise, the expansion-parameter of the Ught-cone correlation function is uM^ rather than 

I^g defined by {u)M^s 



Because of this, in the hght-cone sum-rules we use an "effective" Borel parameter Al^g defined by {u)M^g = 



Mgpj, M2pt being the Borel parameter used in the QCD sum-rules for Jb- 
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rrib [GeV] 


4.7 4.8 4.9 


So [GeV^] 
fB [MeV] 


34.5 ±0.5 33.5 ±0.5 32.5 ± 0.5 

177 ±5 150 ±5 123 ±5 


So [GeV^] 
/b. [MeV] 


35.5 ±0.5 34.5 ±0.5 33.5 ± 0.5 
191 ±5 162 ± 5 135 ± 5 



Table 1: Values for Jb and fs^ from QCD sum- rules in dependent on the 6-quark mass. The 
Borel parameter window is = (4-8) GeV^. 
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Figure 3: Light-cone sum-rule results for 3^ ^ vector-meson form factors. Renormalization 
scale for Ti is = rrih = 4.8 GeV. Further parameters: rrif, = 4.8 GeV, Sq = 33.5 GeV^, 
M2 = 6 GeV2. 
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Figure 4: Light-cone sum-rule results for Bg — > vector-meson form factors. Renormalization 
scale for Ti is fi = nib = 4.8 GeV. Further parameters: nib = 4.8 GeV, sq = 34.5 GeV^, 
M2 = 6 GeV^. 

correlation functions and cancels in the ratio that gives the form factors. Although literally 
we only calculated radiative corrections to the leading-twist contribution to the light-cone 
expansion, it is unlikely that yet unknown corrections to the higher-twist terms could change 
this pattern dramatically. We thus believe that radiative corrections are under good control. 

The next question concerns the convergence of the light-cone expansion. The higher- 
twist terms have several sources: some depend on the intrinsic properties of the multi- 
particle Fock-states of the vector-meson and some appear as meson-mass corrections to the 
two-particle valence state. The latter ones, described in terms of the same parameters as 
the leading- twist distribution amplitudes, turn out to be numerically dominant; this is very 
welcome, as the matrix elements describing the multi-particle states are only poorly known. 
To be specific, putting all intrinsic higher-twist parameters ( in Table |B| to zero, the form 
factors change by at most 3%. Hence, we conclude that the light-cone expansion is under 
good control as well. 
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0.49 


-0.241 
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TP 


0.202 


1.06 


-0.076 


0.261 


1.20 


0.098 


Tf 



Table 2: B^.d-decay form factors in a three-parameter fit. Renormalization scale for Tj 
11 = rrih = 4.8 GeV. The theoretical uncertainty is estimated as 15%. 
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Table 3: iJ^-decay form factors in a three-parameter fit. Renormalization scale for Tj 
ji — rrib — 4.8 GeV. The theoretical uncertainty is estimated as 15%. 
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Figure 5: Separate contributions to the form factors Al and Af-* . 



The dependence of form factors on the sum-rule parameters is small, too. Changing rrib 
by ±100 MeV makes a 5% effect at most and is most pronounced at large g^; at = 
it is a 0.8% effect. This result means that, as for radiative corrections, there is a strong 
cancellation of the rub dependence in the ratio of the light-cone correlation function to fs- 
The same statement holds for the dependence on the continuum-threshold within the limits 
specified in Table |l]. For the dependence on the Borel parameter, we find an ~ 7% effect, 
increasing with g^, which again reminds us of the fact that the light-cone sum rules become 
less reliable for large g^. 

The overall normalization of the form factors depends on the vector-meson decay con- 
stants fv and fy, the former determined experimentally, the latter calculated from QCD 
sum rules (see Table ^). The corresponding uncertainty is at most 3%. 

Adding up all the errors in quadrature, we obtain an uncertainty of the form factors of 
-11%. 

The shape of leading-twist distribution amplitudes, characterized by the Gegenbauer 
moments af^ for the p and af'^., cl^'k* for the K*, affects most significantly the slope of 
the form factors and is illustrated in Fig. ^ by two examples: and A^* . The curves 
labelled "asymptotic" designate the form factors as obtained by putting the to zero in 
Eqs. ( p.l5|) and ( p.30|) ; the corresponding meson distribution amplitudes are completely 
model-independent and dictated by perturbative QCD. The curves labelled "oj" show cor- 
rections to this limit, which take into account non-perturbative corrections to the distribution 
amplitudes; for illustration we assumed in this figure the value a 



4,P 



4,p 



0.1 at /i = IGeV 

as a ball-park estimate for potential higher-order terms; this contribution is not included in 
the final results. The curves labelled "twist-3 and 4 terms" show the contribution induced 
by the C's in Table ^ and for the K* also contain terms explicitly proportional to the strange 
quark-mass. It is obvious that the "asymptotic" contribution grossly dominates, and the re- 
maining terms only add marginal corrections. It is also obvious that the twist-3 and 4 terms 
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F 


Ai 


A2 


V 


Ti 


T3 


F^*(0)/F''(0) 


1.30 ±0.13 


1.28 ±0.13 


1.36 ±0.14 


1.33 ±0.13 


1.29 ±0.13 



Table 4: Size of SU(3)-breaking for i?„_rf-decays into p or K*. 



do not have much overall influence, whereas the contribution in 02 (for A^) and ai (for Af ) 
tend to slow down the increase of the form factors as functions of q^. All involved parameters 
(except the couplings fy and fy) come with considerable theoretical uncertainty. However, 
the only important error is that in a2,p and ai^K*'- it contributes ~ 10% to the uncertainty in 
our predictions. Adding this number (in quadrature) to the ~ 11% error from other sources, 
we end up with a total uncertainty of light-cone sum-rule predictions of order ~ 15%, which 
is our final error estimate. An improvement is to be expected if future lattice calculations 
achieve an accuracy better than that quoted in Table 0. 

A few remarks are in order on the pattern of SU(3)-breaking. It enters our calculation 
at the following places: 

• difference in decay constants: fx^/fp ~ fK*/fJ = 1-14, IbJIb = 1-08; 

• different meson-masses and continuum-thresholds sq (Table 0); 

• different vector-meson distribution amplitudes (Table 

Figure ^ also illustrates the relative size of these effects: the difference between the "asymp- 
totic" curves is almost exclusively due to the difference in fp and Jk* and makes a 17% 
effect. For K*, the 02 are small, whereas the ai are large and thus increase the form factor. 
For Bs K* decays, the sign in ai is negative and Jb, is larger than Jb, so that we observe 
considerably smaller form factors, see Table |^. The total SU(3)-breaking corrections amount 
to ~ 35%, half of which comes from the decay constants and half from the bigger momentum 
carried by the s quark in the strange hadron. Specifically, for Bu^d-decaj form factors at 
= we obtain the values given in Table ^. 

In Fig. 1^ we present a comparison of our results for B p semileptonic and rare radia- 
tive form factors with the lattice calculations by the UKQCD collaboration The 
agreement is very good. We wish to emphasize that the light-cone sum-rule approach is 
theoretically more sound at small values of g^, and in this sense is complementary to lattice 
techniques, which work best in the large region. A similar comparison for B —>■ K* decays 
is presented in Fig. |^. The agreement is somewhat worse in this case; the lattice data favour 
smaller SU(3)-breaking effects. This question deserves further study. 

Finally, in Table |^ we present a comparison of the results of this work for the form factor 
values at g^ = with earlier sum-rule calculations and the lattice results obtained using the 
light-cone sum-rule constraints. 
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^2[GeV^] 

Figure 6: Comparison of the light-cone sum-rule predictions for the B p form factors with 
lattice calculations ||22|, |2^. Lattice errors are statistical only. The dashed curves show the 
15% uncertainty range. 
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Figure 7: Comparison of the light-cone sum-rule predictions for the B K* form factors 



with lattice calculations [^2] , |23[| . Lattice errors are statistical only. The dashed curves show 
the 15% uncertainty range. 
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Table 5: Comparison of results from different works on form factors at = 0. 

4 The Heavy- Quark Limit 



The behaviour of 5-decay form factors in the limit rub 



oo is interesting for various reasons. 

so that in this paper we 



This limit was already discussed in some detail in Refs. ||18| , |T2, 
only summarize the main points. 

The first question concerns the scaling behaviour of form factors as functions of the b- 
quark mass. The behaviour depends on the momentum-transfer and is different for small and 
large recoil. For 0, or, more precisely, for ml — ~ O(m^), all form factors in question 

3 /2 

scale as ~ ^/^b ■ This behaviour can be proved in perturbative QCD, taking into account 



Sudakov suppression of large transverse distances, but it is not restricted to this regime 
and extends to "soft" terms as well [^, For ml — ~ 0{mh), on the other hand, the 
form factors obtained from light-cone sum-rules satisfy the scaling behaviour predicted by 
Heavy-Quark Effective Theory (HQET) ^5[. For realistic values of the 6-quark mass, these 
two regimes are not well separated; therefore large corrections to asymptotic scaling are to 
be expected. Some estimates of pre-asymptotic corrections are presented in Refs. |TB|, |I^]. 
They have to be considered as indicative only. We do not attempt to further quantify 
pre-asymptotic corrections in this work. 

The second question concerns possible relations between different form factors in the 
heavy-quark limit. Heavy-quark symmetry implies exact relations between semileptonic and 
penguin form factors at small recoil and renormalization scale fi = m^ p5|, which, using the 
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penguin form factor definitions in Eq. ( |2.13| ), can conveniently be written as: 



A{q') - B{q') 
C{q') 



2m 



B 



ms + my 



V{q' 



{ruB + mv) ^^^^2^ _ m^B - q'^ + ^(?^) 



niB 



ms - my y^^-^ ^ vrLy(m\ - m\) 



rriB 



^B - my 2 

H M{q 



TUB niB+my 

[Aoiq') - Asiq')] 



(4.1) 
(4.2) 



TUBq 



2m 



B 



(4.3) 



Writing tlie relations in this form emphasizes their different behaviour in the heavy-quark 
limit: at small recoil, both sides of Eqs. (|4.1| ) and (|4.3| ) are of order A/mt, while Eq. (|4.2| ) 
relates combinations of form factors, which are of order 1 / ^/m^,. The numerical comparison 
for B ^ p transitions is presented in Fig. ^. We note that (a) Eq. ( |4.1| ) is satisfied with 



high accuracy, (b) The relation (4^) is violated. However, both sides are numerically small 
compared to Eq. ( |4.1| ), in agreement with the expected 1/mf, suppression, (c) The relation 
( [4.3|) is very well satisfied at — >• and it holds with 20% accuracy at large g^; both sides 
turn out to be small at large recoil, which implies significant cancellations between the terms 
on the right-hand side. 

For phenomenological applications it is more appropriate to rewrite the Isgur-Wise rela- 
tions (|4.1| )-( |^I^ ) in terms of the form factors defined in ( |1.3[ ): 



- my 
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(4.4) 



(4.5) 



+ q^ 



m 



V 



2m Bq^ 



{mB + my)Ai{q^) - (m^ - my)A2{q^) .(4.6) 



Note that such a rewriting mixes terms of different order in 1/ m^, in the small recoil region, 
and in this sense is not fully consistent with the derivation in pSl. It can be justified. 



however, by observing that the hierarchy of contributions is different at large recoil and 
all the terms become formally of the same order. The numerical comparison for B ^ p 
transitions is presented in Fig. |^. The accuracy proves to be excellent for the relation ([4.4|) , 
which is observed to within 3% accuracy, and good for ( |4.5D with deviations of at most 8%. 
Relation (|4.6|), however, is violated by 20% for q^ > 15 GeV^. Since fidelity of the sum-rules 
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Figure 8: Isgur-Wise relations ( |4.1| )-( P7^ ) for B ^ p transitions. Renormalization scale is 
p = mi,. Solid and dashed curves correspond to expressions appearing on the l.h.s. and r.h.s., 
respectively. 
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Figure 9: Isgur-Wise relations ( |4.4| )-( |4^ ) ioi B ^ p transitions. Renormalization scale is 
fi = mi,. Solid and dashed curves correspond to expressions appearing on the l.h.s. and r.h.s., 
respectively. 
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worsens in the high-g^ region, it is not clear whether this disagreement indicates a genuine 
1/mf, correction or is an artefact. Our results reinforce an earlier observation in [|l^ that 
the relation in ( |4.4| ) is satisfied within ~ (5-7)% in the whole region of to leading-twist 
accuracy in the light-cone sum-rule approach, and strongly support the conjecture of p6 
about the validity of heavy-quark symmetry relations in the region of small in heavy-to- 
light transitions. 



5 Conclusions 

We have given a complete analysis of 5-decay form factors to light vector-mesons in the 
light-cone sum-rule approach. The principal new contribution of this work is the calculation 
of radiative corrections and higher-twist corrections to the sum-rules, which are calculated 
for the first time. We observe that the light-cone sum-rules turn out to be very robust against 
corrections in the light-cone expansion, whose numerical impact proves to be minimal. Ra- 
diative corrections seem to be well under control. In cases where higher-twist corrections are 
important, they are dominated by meson- mass effects, which do not involve free parameters. 
The theoretical accuracy of the approach is thus restricted entirely by the duality approx- 
imation for the extraction of the i?-meson state from the continuum and contributions of 
higher resonances. The usual "educated guess" is that the accuracy of such an extraction is 
of order 10%, which provides an irreducible error. Effects of yet higher radiative corrections 
and yet higher-twists are likely to be much less; therefore, the sum-rules derived in this 
work cannot be improved significantly. The numerical analysis, however, can and should 
eventually be updated, once estimates for the meson distribution amplitudes, 6-quark mass 
and fB become more precise. In particular, lattice calculations of the tensor couplings fy 
and the parameters a}\ 21 ^12 meson distribution amplitudes would be most welcome. 
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A One-loop Integrals 

For the calculation of radiative corrections, 

(27r)^ {k + upYk^[{q - kf - m?] " 

[ d^k ka{q - k)p _ 

J [27^)'^ {k + upfk^[{q - kf - m^] ~ 



we need the following integrals: 
Aupo, + Bq^, (A.l) 

Cqc^p + Dq^qp + Eq^upp + Fup^qp + . . • ,(A.2) 
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(27r)-° k'^{k — upY{k + upY[{up -\- q — ky — w? 
d^k 1 



{2t:)^ k'^[{up + q - kY - m?] 
{2t:)^ k^[{up + q - kf - m?\ 



^ Z {q + up)c 



Hgafs + Iqaqp + (A.3) 

(A.4) 
(A.5) 



where the dots stand for terms that are irrelevant to the present calculation. The functions 
A, B, C are obtained from A, B and C by the replacement 



u ^ u\ 



-PB- 



We shall use the notation 



s = rr? — upg — uf^ ^ - = 7e + In Att, 

e e 



(A.6) 



(A.7) 



with D = 4 — 2e. In order to perform Borel transformation and continuum subtraction, the 
following spectral representations for the above integrals prove the most convenient: 
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B Summary of Meson Distribution Amplitudes 



The expressions collected in this appendix are principally the result of recent studies reported 
in Refs. 0, |10|, 0]. We use a simplified version of the set of twist-4 distributions |]ll|, taking 
into account only contributions of the lowest conformal partial-waves, and for consistency 
discard contributions of higher partial- waves in twist-3 distributions in cases where they enter 
physical amplitudes multiplied by additional powers of mp. The SU(3)-breaking effects are 
taken into account in leading- twist distributions and partially for twist-3, but neglected 
for twist-4. Explicit expressions are given below for a (charged) p-meson. Distribution 
amplitudes for other vector-mesons are obtained by trivial substitutions. 

Throughout this appendix we denote the meson momentum by and introduce the 
light-like vectors p and z such that 



(B.l) 



The meson polarization vector e^^^ is decomposed in projections onto the two light-like 
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vectors and the orthogonal plane as 

(eW ■ z) 



,(A) 



pz 



m 



2pz 



(B.2) 



B.l Chiral-even distributions 

Two-particle quark-antiquark distribution amplitudes are defined as matrix elements of non- 



local operators on the light-cone |]10 

{Q\u{z)^^d{-z)\p-{P,\)) = f,mp 



P-z Jo 



2^{p-zY 



1 e'-'**^ ■ z 

-z„- —mi I due'^^ '^gs^u, f/) 
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and 



(0|«(.)7,75C^(-.)|p-(P, A)) = i f/p - /j!!^!i±^ 



rrir 



Jo 



(B.4) 

For brevity, here and below we do not show the gauge factors between the quark and the 
antiquark fields and use the short-hand notation 

^ = u-{l-u) = 2u-l. 
The vector and tensor decay constants fp and fj are defined, as usual, as 

(0|^(0)7;.rf(0)|p-(P,A)) = /pm.eW, (B.5) 

{0\u{0)a,,d{0)\p'iP,\)) = z/J(eWp, -eWp^). (B.6) 

The distribution amplitude is of twist-2, g^^^ and g^^^ are twist-3 and g^ is twist-4. All 
four functions = {0||, g^^\ g^_^\ gs} are normalized as 



du(j){u) = 1, 



(B.7) 



which can be checked by comparing the two sides of the defining equations in the limit 
Zf^ ^ and using the equations of motion. We keep the (tiny) corrections proportional to 
the u and d quark-masses and rrid to indicate the SU(3)-breaking corrections for K*- and 
0-mesons. 

In addition, we have to define three-particle distributions: 



{Q\u{z)gG^^-fa75d{-z)\p (P,A)) = fpmpPa\pye^^l 



Pi^e^il]^{v,pz) 
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,(A) . 



+ fpmp ^^^^^ Pa[Pf,Zu - PuZ^,]'^{v,pz) (B.8) 



where 



{0\u{z)gGf,^i-fad{-z)\p (P)) = fpmppj^p^e^^l-p^e^llv{v,pz) 

. e(^) ■ z 
pz 

.e(-^) • z 

+ fpmp-^—^Pa\p^,z^-VuZ^]^{v,pz), (B.9) 

A{v,pz) = j Pae-^P"(""-"'*+^"«U(a), (B.IO) 

etc., and a is the set of three momentum fractions a = {ad, ctu, C(g}- The integration measure 
is defined as i i i 

J Va = J dad J dau J dag 6{1 — aj). (B-H) 

The distribution amphtudes V and A are of twist-3, while the rest is twist-4 and we have 
not shown further Lorentz structures corresponding to twist-5 contributions^. 

Calculation of exclusive amplitudes involving a large momentum-transfer reduces to eval- 
uation of meson-to-vacuum transition matrix elements of non-local operators, which can be 
expanded in powers of the deviation from the light-cone (see text). To twist-4 accuracy one 
can use the expression for the axial- vector matrix element in ( |B.4| ) as it stands, replacing the 
light-cone vector by the actual quark-antiquark separation x^. For the vector operator, 
the light-cone expansion to the twist-4 accuracy reads: 



I e^'^^x r rn?x'^ 

(0|M(x)7^rf(-x)|p-(P, A)) = /pmJ -— / due'^^"" (p\\{u, IJ.) + — ^A(M,/i) 







where 



rfMe^«-f'^C(n,p)|, (B.12) 
C(n) = g^{u) + <P\\{u) - 2g^l\u) (B.13) 



e^^^x 2 



^ We use a different normalization of three-particle twist-3 distributions compared to pOl 
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and A(m) can be related to integrals of three-particle distributions using equations of motion. 
All distribution functions in (p.l2| ) are assumed to be normalized at the scale //^ ~ (to 
leading- logarithmic accuracy). In practical calculations it is sometimes convenient to use 
integrated distributions 



C»(m) = - / dvC{v), = - / dvC'^'^v) 



(B.14) 



For the leading twist-2 distribution amplitude 0|| we use 

3 



(j)\\{u) = 6uu 



l + 3a\^ + al^{5e-r 



[B.15) 



with parameter values as specified in Table ^ The expressions for higher-twist distributions 
given below correspond to the simplest self-consistent approximation that satisfies the QCD 
equations of motion [|10], pi]] : 



Three-particle distributions of twist-3: 



A{a) = 360 Csadttual 1 + uf^Tag - 3)] 



[B.16) 
[B.17) 



Two-particle distributions of twist-3: 

g^f\u) = Quu 



-1-6 5+ {3uu -|- M In M -|- M In m) -|- 6 5- (m In M — -u In m) 



(B.18) 



+ 



112 



(3 - 30e^ + 35^^) 



3 ~ 3 ~ 

+ -6+{2 + \nu + lnu) + -5^{2^ + lnu- In u) , 



(B.19) 



Three-particle distributions of twist-4: 



30(1 — ag)a 



$(a) 



- 3C3 + 3(4^ 



30(a„ — ad)a 
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= 0. 



2 1 ■ 

3C3 + 3C4 



(B.20) 



Two-particle distributions of twist-4 

A{u) -- 



4 20 8 

5 + Y^^ + 9^' 



g^{u) = 6m(1 -u) + 

C{u) -- 



10 20/ 
y C4 - yCs 



:i-3a, 



3 10 10 

2 + y + y 



;i - 3r 



3 10 10/ 

2 + y + y 



(B.21) 



where the dimensionless couphngs Cs and (4 are defined as local matrix elements 



{0\ugGp^-faj5d\p {P, X)) = fpUipCs 



1 



+ 3/p"^pC4 



(B.22) 



and have been estimated from QCD sum-rules |23, |2^. The terms in 5± and 5± specify 
quark-mass corrections in twist-3 distributions induced by the equations of motion. The 
numerical values of these and other coefficients are listed in Tables and 0]. Note that 
we neglect SU(3)-breaking effects in twist-4 distributions and in gluonic parts of twist-3 
distributions. 

B.2 Chiral-odd distributions 

There exist four different two-particle chiral-odd distributions defined as 



{0\u{z)ap,d{-z)\p-{P,X)) = ifl 



e(^) ■ z 



(p-/2 P 



ml / due'^P-'h'^\u,p^) 



P z Jo 

^In the notation of Ref. [[0|, wJ^q = '^3^. C3 = Cs, and (Y = (3/28)C3W3' 



due'^P-'h:i{u,fi'^ 



(B.23) 
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{Q\u{z)d{-z)\p-{P,X)) = -z (fj - f^Hhil^) (e(A) . ,)ml f due^^^^hf {u, f,'). (B.24) 

\ "^P J Jo 

The distribution amplitude (f)± is twist-2, h\^^^ and /ly*^ are twist-3 and is twist-4. All four 

functions (f) = h^^\ h^i^\ h^} are normalized to j^du(j)[u) — 1. 

Three-particle chiral-odd distributions are defined to twist-4 accuracy as 

{^\u{z)a^figG^^y{vz)d{-z)\p~{P,\)) = 

e(^) ■ z 

, rT 2\ (A) ± (A) ± (A) ± , (A) ± nn(4)/ n 



+ 



p p\ 



\PaP,xe^±0Z^ - P^i.e'llz^ - PaPye'-llz^,+ ppp^e'-llz^]T^^^\v,pz) 



(A)^ 



+ 



pz 
pz 



IPaP^e^^^Up - ppPi^e^^lza - PaPue^±!,Zfs + ppp^e^^l^Za]T^^' {v,pz) 



(A), 



,(A)^ lm(4). 



+ 



(B.25) 



and 



{Q\u{z)gG^,{vz)d{-z)\p-{P, A)) = if,ml[ef^p, - e^^ip,]S{v,pz), 

{Q\u{z)igG,,{vz)^,d{-z)\p-{P, A)) = ifjmy^y - e^^lp,]S{v,pz). (B.26) 

Of these seven amplitudes, T is twist-3 and the other six are twist-4. 

The light-cone expansion of the non-local tensor operator can be written to twist-4 ac- 
curacy as 



{0\u{x)a^,d{-x)\p-{P,X)) = 

-ifjliei^^P^-el'^P,) fdue^^^^ 
L -'0 



2 2 
TTl X 

0±W + ^AT(ti) 



= (A) . 



X 



{Pxf " 

2 ,.1 



du(^^^''^T{u) 



(B.27) 
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where and Cr are expressed in terms of the distribution amphtudes defined above as 



i>t{u) = h^^\u) - - ^h{u), 



Ct{u) = hs{u) - 



(B.28) 



and At can be related to integrals of three-particle distribution functions using the equations 
of motion. 

We introduce the notation, similar to Eq. ([B.14|): 



u 



dvMT{v), 



dv Ct{v). 



For the leading twist-2 distribution amphtude (j)± we use 

3 



(f)±{u) = Quu 



l + 3ai^ + ai-{5e~l) 



(B.29) 



(B.30) 



with parameter values as specified in Table ^ The expressions for higher-twist distributions 
given below correspond to the simplest self-consistent approximation that satisfies all QCD 
equations of motion [|10], |ll[|: 

• Three-particle distribution of twist-3: 



T{a) = 540C3uJ{ad - au)adaual. 



(B.31) 



Two-particle distributions of twist-3: 



h^^^\u) = Quu 



-|- 3 (5+ (3mm -|- m In m + m In m) + 3 (5_ (m In u — m In u) , 



(B.32) 



hf{u) = ^e+l ai - 1) + e {5e - s) + ^C3^3^(3 - 3oe + 35^) 



3 3 
+ - 5+ (1 + ^ \nu/u) + - (5„ ^ (2 + Inw + Inw) 

Three-particle distributions of twist-4: 

Ti'\a)=Tt\a) = 0, 
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(B.33) 



V 




K* 


K* 





J V [-LVJ-*-- V J 


198 ± 7 


226 ± 28 


226 ± 28 


254 ± 3 


/f[MeV] 


160 ± 10 
152 ± 9 


185 ± 10 
175 ± 9 


185 ± 10 
175 ± 9 


215 ± 15 
204 ± 14 


a! 





0.19 ±0.05 
17 ± 04 


-0.19 ±0.05 
-0 17± 04 





a| 


0.18 ±0.10 
16 ± 09 


0.06 ±0.06 
05 ± 05 


0.06 ±0.06 
05 ± 05 


0±0.1 







0.20 ±0.05 
0.18 ± 0.05 


-0.20 ±0.05 
-0.18 ± 0.05 







0.20 ±0.10 
0.17 ± 0.09 


0.04 ±0.04 
0.03 ± 0.03 


0.04 ±0.04 
0.03 ± 0.03 


0±0.1 







0.24 
0.22 


0.24 
0.22 


0.46 
0.41 







-0.24 
-0.22 


0.24 
0.22 





s+ 





0.16 
0.13 


0.16 
0.13 


0.33 
0.27 


5_ 





-0.16 
-0.13 


0.16 
0.13 






Table A: Masses and couplings of vector-meson distribution amplitudes, including SU(3)- 
breaking. In cases where two values are given, the upper one corresponds to the scale //^ = 
1 GeV^ and the lower one to = m| — ml = 4.8 GeV^, respectively. We use 771.5(1 GeV) = 
150 MeV and put the u and d quark mass to zero. 





Cs 






T 


C4 


CI 


CI 




0.032 


-2.1 


3.8 


7.0 


0.15 


0.10 


-0.10 


V 


















0.023 


-1.8 


3.7 


7.5 


0.13 


0.07 


-0.07 



Table B: Couplings for twist-3 and 4 distribution amplitudes for which we do not include 
SU(3)-breaking. Renormalization scale as in the previous table. 
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S{a) = 30Cj(l - ag)al, 



S{a) = 30Cj(l - ag)al. (B.34) 



Two-particle distributions of twist-4: 

hiu) = 6n(l -u) + 5[cr + cr](i - 3e), 

^2 4 



At{u) = 30u\l - uy 



5 + 3^^^ 3^^"^ 



[B.35) 



The constants (J and (J are defined as 

{0\ugG,,d\p-{P, A)) = 2/Jm3cJ(eWp, - e^P^), 

(0|M^G^,^75C^|p-(P, A)) = tfJmlClief^Pu - e^^^P^) (B.36) 

and liave been estimated in |]l| from QCD sum-rules: 

Cr~-Cj~0.10. (B.37) 

Other parameters are given in Table [Afl . As in the chiral-even case, we neglect SU(3)- 
breaking corrections in twist-4 distributions. 
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